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^ Abstract 

Low-energy baryon-baryon interactions are calculated in a high-statistics lattice QCD study on a 
single ensemble of anisotropic clover gauge-field configurations at a pion mass of m n ~ 390 MeV, a 
spatial volume of L 3 ~ (2.5 fm) 3 , and a spatial lattice spacing of b ~ 0.123 fm. Liischer's method 
is used to extract nucleon-nucleon, hyperon-nucleon and hyperon-hyperon scattering phase shifts 
at one momentum from the one- and two-baryon ground-state energies in the lattice volume. The 
isospin-3/2 AS interactions are found to be highly spin-dependent, and the interaction in the 3 Si 
channel is found to be strong. In contrast, the iVA interactions are found to be spin-independent, 
within the uncertainties of the calculation, consistent with the absence of one-pion-exchange. The 
only channel for which a negative energy-shift is found is AA, indicating that the AA interaction 
is attractive, as anticipated from model-dependent discussions regarding the H-dibaryon. The NN 
scattering lengths are found to be small, clearly indicating the absence of any fine-tuning in the 
NN-sector at this pion mass. This is consistent with our previous Lattice QCD calculation of NN 
interactions. The behavior of the signal-to- noise ratio in the baryon-baryon correlation functions, 
and in the ratio of correlation functions that yields the ground-state energy splitting is explored. 
In particular, focus is placed on the window of time slices for which the signal-to-noise ratio does 
not degrade exponentially, as this provides the opportunity to extract quantitative information 
about multi-baryon systems. 
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I. INTRODUCTION 



The strong interactions among baryons are key to every aspect of our existence. The two- 
and higher-body interactions among protons and neutrons conspire to produce the spectrum 
of nuclei and the complicated chains of nuclear reactions that allow for the production of 
the elements forming the periodic table at the earliest times of our universe, in the stellar 
environments that follow, and in reactors and our laboratories. Decades of experimental 
effort have provided a very-precise set of measurements of the nucleon-nucleon scattering 
cross sections over a wide range of energies pQ, and these cross sections have in turn given rise 
to the modern nuclear forces. These experimentally determined two-body forces, encoded 
by potentials such as AVi 8 [2] and the chiral potentials [3], when supplemented with three- 
body interactions, now provide the cornerstone of our theoretical description of nuclei [I] and 
their interactions. The two-nucleon forces are observed to be significantly more important 
than the three-nucleon forces, which, in turn, are significantly more important than the 
four-nucleon forces. Present-day calculations are sufficiently precise that the inclusion of 
three-nucleon forces is required in order to post-diet the structure of light nuclei 016]. Given 
the relatively small contributions of the three-nucleon and higher-body interactions to light- 
nuclei, where reliable calculations are presently possible, there is considerable uncertainty 
with regard to their form. This leads to enhanced uncertainties in the calculation of systems 
for which there is little or no experimental guidance, such as moderate to high-density 
neutron-rich environments, and more generally, nuclear environments at densities exceeding 
that of nuclear matter. 

In dense nuclear systems it is not only the multi-nucleon forces that are difficult to 
quantify as non-nucleonic objects may play an important role. In a core-collapse supernova, 
it is the nuclear equation of state (NEOS) that ultimately dictates whether the system 
collapses into a neutron star or forms a black hole. This, in turn, is determined by the 
composition and structure of the hadronic matter in the core of the supernova, which is at 
(baryon number) densities that are a few times that of nuclear matter. At such densities, the 
strange quark may play a pivotal role through the formation of a charged kaon condensate 
made possible by the strength of attractive kaon-nucleon interactions. It may also become 
energetically favorable for the matter to contain strange baryons, such as the £~ or A, 
because of their interactions with nucleons. In either case, it is the two-body interactions 
between the strange hadrons and the nucleons, and also between themselves, that largely 
determine the composition of the hadronic matter at core-collapse densities, and ultimately 
the fate of the collapsing supernova [7] . Unfortunately, experimental determinations of the 
interactions of strange hadrons are very challenging due to their weak decays, and existing 
cross-section measurements are not precise enough to provide meaningful constraints on the 
NEOS at core-collapse densities. So, while it is the interactions between three-nucleons that 
currently limit the precision with which properties of material composed of neutrons and 
protons can be calculated, it is the two-body interactions involving strange hadrons that 
currently provide the most serious limitation to reliable calculations at densities that exceed 
that of nuclear matter. 

One of the major objectives of lattice QCD (LQCD) is to calculate the properties and 
interactions of nucleons and, more generally, systems comprised of multiple hadrons. Precise 
exploration of the simplest multi-hadron systems has recently become possible with signif- 
icant advances in computational resources, as well as through algorithmic and theoretical 
developments (for a recent review, see Ref jS]). The 1 = 2 two-pion system, 7r + 7r + , is the 
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simplest of such multi-hadron systems to calculate in LQCD, and current computational re- 
sources have allowed for a ~ 1% level determination of the 7r + 7r + scattering length [9| [T0llTT] . 
Further, systems comprised of up to twelve 7r + 's [121 EE] and systems comprised of up to 
twelve K +, s [H] have been explored, allowing a determination of the three- ir + and three- K + 
interactions as well as aspects of pion and kaon condensates. 

In general, a determination of the two-particle scattering amplitude, or multi-body in- 
teractions, with LQCD requires calculating the energy eigenvalues of the appropriate two- 
hadron system in the finite volume [151 EH HZ]- The energy differences between the multi- 
particle energy levels in the finite volume and the sum of the particle masses determines 
the scattering amplitude at the corresponding energy. Processes of interest to low-energy 
nuclear physics occur in the MeV energy regime, while the masses of the baryons and nuclei 
are in the GeV regime. As a result, very precise (high-statistics) measurements must be per- 
formed in order to reliably extract useful constraints on two- and higher-body interactions 
from a lattice calculation. 

In contrast to mesonic systems, lattice QCD calculations in the baryon sector are com- 
plicated by the statistical behavior of the correlation functions |18j . We performed the first 
unquenched LQCD calculation of nucleon-nucleon (NN) [19] and hyperon-nucleon (YN) [20] 
scattering lengths, and found that the NN scattering lengths appear to be of natural size at 
larger pion masses, and are set by the range of the interaction. Consequently, an increase in 
the magnitude of the scattering length as the pion mass is reduced down to its physical value 
is anticipated. Further, we have recently performed the first lattice QCD calculation of scat- 
tering in various meson-baryon channels [21] using the mixed-action scheme of domain-wall 
valence quarks on the MILC staggered sea. 

Recently, we have performed a high-statistics calculation of a number of single-hadron 
correlation functions [22] and also the first calculations of three-baryon systems [23] on an 
ensemble of the anisotropic gauge-field configurations generated by the Hadron Spectrum 
Collaboration [24] 125] . This ensemble has a pion mass of m n ~ 390 MeV, a spatial lattice 
spacing of b s = 0.1227 ± 0.0008 fm, an anisotropy £ = b s /b t = 3.500 ± 0.032 and a lattice 
volume of L 3 x T = 20 3 x 128. The goal of those studies was to "jump" an order of magnitude 
in the number of measurements performed to estimate correlation functions, and to explore 
the "new territory" that subsequently emerged. The baryon masses were extracted with fully 
quantified uncertainties at the < 0.2%-level from the 0.28 x 10 6 measurements performed 
on 1194 gauge-field configurations. With a somewhat smaller statistical sample, the binding 
energies of the pnn (triton = 3 He by the isospin symmetry of the LQCD calculation) and 
"H°S°n" three-baryon systems were investigated and in the latter case, determined with an 
uncertainty of ~ 3 MeV/baryon. Subsequently, a study of the 3 He (ppn) and 4 He (ppnn) 
systems in quenched QCD at a large quark mass appeared [26]. A number of important and 
surprising observations were made in those works [221 1231 12S] that have modified how we 
foresee moving toward calculating the properties of light nuclei. 

One of the most important aspects of our previous work was the detailed study of the 
signal-to-noise ratio in the single-baryon correlation functions. The signal-to-noise ratio was 
found to be approximately independent of time for a significant number of time-slices [221 E3] 
prior to evolving toward the expected exponential degradation [18J. This window of time- 
slices is understood in terms of the suppression of contributions from purely mesonic states 
in the correlation function that determines the variance of the single-baryon correlation 
function. Given that the signal-to-noise ratio for a system containing more than one baryon 
is expected to scale (approximately) as the product of the signal-to-noise ratio's of the 
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individual baryons (neglecting their interactions), this window of time slices suggests that 
calculations of the energy levels of systems containing a number of baryons in this lattice 
volume with these interpolating operators is possible. In this work we present the baryon- 
baryon scattering phase shifts that have been measured in our high-statistics anisotropic 
clover-quark calculation. 

Our analysis and results are presented in the following manner. In Section [XXJ, we introduce 
the details of the formalism used in extracting phase shifts from our LQCD calculations. 
Section III briefly reviews our single baryon results before we present our extractions of 
two-baryon interactions in Section IV Section [V] details our analysis of statistical scaling 
and noise in the various correlation functions while Section [VI] is a concluding discussion. 



II. LATTICE QCD CALCULATIONS 

A. Liischer method for extracting scattering parameters 

In this work, the finite volume scaling method (Luscher's method) [151 EE E] is employed to 
extract the two-particle scattering amplitudes below inelastic thresholds at a given energy. 
In the situation where only a single scattering channel is kinematically allowed, the deviation 
of the energy eigenvalues of the two-hadron system in the lattice volume from the sum of 
the single-hadron masses is related to the scattering phase shift, 5. For energy eigenvalues 
above kinematic thresholds where multiple channels contribute, a coupled-channels analysis 
is required as a single phase shift does not parameterize the S-matrix. The energy shift for 
two particles A and B, AE = Eab — Ea — Eb, can be determined from the correlation func- 
tions for systems containing one and two hadrons. For baryon-baryon systems, correlation 
functions of the form 

C B;r (p,t) = Yl ^ r « ^«( x '*) ^( x o>°)> (!) 

X 

C BllB3;r (Pi,P2,f) = e ipi - xl e iM ^r^(B liai (x 1 ,OB2,« 3 (x 2 ,t)Bi A (xo,0)H 2A (x ,0)), 

XI, X2 

are used, where B denotes a baryon interpolating operator, and (3^) are Dirac indices, 
and the Y are spin tensors that typically project onto particular parity and/or angular 
momentum states. The baryon octet interpolating operators are of the form 

p«(x,t) = e^<(x,t)(^ T (x,t)C 7 5^(x,t)) , 

A a (x,t) = e^V a (x J t)(^" r (x,t)C r T5 d*(x,t)) , 

E+(x,t) = e^<(x,t)(^ T (x,t)C 7 55 fc (x,t)) , 

S°(x,t) = e^V a (x,t) (^ T (x,t)C 75 s fc (x,t)) , (2) 

where C is the charge-conjugation matrix and ijk are color indices. Other hadrons in 
the lowest-lying octet can be obtained from the appropriate combinations of quark fla- 
vors. The parentheses in the interpolating operators indicate contraction of spin indices 
into a spin-0 "diquark". It is worth pointing out that the overlap of the composite oper- 
ator jBi )Q1 (xi,t)^2,a 2 ( x 25^) on t° two-baryon states is not simply the product of individual 
baryon "Z-factors", but depends explicitly upon xi — x 2 , with a correlation length set by the 
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pion mass. This fact precludes a determination of interpolator-independent baryon-baryon 
potentials from lattice QCD calculations |27j . 

Away from the time slice on which the source is placed (in this case t = 0) these correlation 
functions behave as 

Cn A (p,t) = ^a(p) Z$(P) e-^W (3) 

n 

C HAHB ( P ,-p,t) = £ Z%ab(p) ^vL(p) e~ EkAB) W \ (4) 



with E { A) {0) = m A and the E {AB) (0) are the energy-eigenvalues of the two particle system 
(we only present calculations of two-baryon systems for which the center-of-mass is at rest) 
in the lattice volume. At large times, the ratio 

C HA (0,t)C HB (0,t) ^ Z oMP)Z OAB ( P ) e (5) 

decays as a single exponential in time with the energy shift, AE (AB) (0). In what follows, 
only the case p = is considered. The energy shifts 



AE (ab) = AE^ B \0) = Ei AB \0) - mA - mB = yj q * + m 2 + ^ q 2 + m | _ m ^ _ mfi 

= 2^- + -. (6) 

(where ^ab = ^ J 4 m s/( m A + m _B) is the reduced mass of the two-particle system) determine 
squared momenta, (which can be either positive or negative). Below inelastic thresholds, 
these are related to the real part of the inverse scattering amplitude via 1 

q n cot 5(q n ) = Lfc (J^J \ , (8) 

where 

|j|<A 1 

S(x) = lim V — 4vr A , (9) 

j IJI 

thereby implicitly determining the value of the phase shift at the energy AE iAB) (or center 
of mass momentum q n ), 8(q n ). Thus, the function pcot 5(p), that determines the low-energy 



1 Calculations performed on anisotropic lattices, such as those used in this work, require a modified energy- 
momentum relation, and as a result eq. ([6]) becomes 

AEi AB ~> ee Ei AB > -m A -m B = Jq*/£ A + m A + y/ql/O+rr^ - m A - m B , (7) 

where are the anisotropy factors for particle A and particle B, respectively, determined from the 
appropriate energy-momentum dispersion relation. The masses and energy splitting are given in terms of 
temporal lattice units and q n is given in spatial lattice units. In the present work we find that the various 
£4 agree within uncertainties and have used £ A = £b = £, for all scattering processes. 
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elastic-scattering cross-section, A(p) oc (p cot 8{p) — ip) 1 , is determined at the energy 
AE^ AB \ 

For a scattering process for which the exchange of a single pion (One-Pion-Exchange 
(OPE)) is allowed by spin and isospin considerations, the function pcot5(p) is a analytic 
function of |p| 2 for |p| < m n /2 (determined by the t-channel cut in the scattering amplitude). 
In this kinematic regime, p cot 5(p) has a series-expansion (the effective range expansion) of 
the form 

pcotS(p) = — + - r |p| 2 + ... , (10) 

where a is the scattering length (with the nuclear physics sign convention) and ro is the 
effective range. While the magnitude of the effective range (and higher terms) is set by the 
pion mass, the scattering length is unconstrained. For scattering processes where OPE is 
not allowed, the lower limit of the cut in the t-channel and the location of inelastic threshold 
set the radius of convergence of the effective-range expansion of pcot5(p). 



B. Expectations at finite temporal extent 

For baryon-baryon correlation functions computed from quark propagators that are anti- 
periodic in the time direction, and using the sources of eq. ([2]), the correlation functions 
contain contributions from hadronic states propagating forwards and backwards in time. 
They are expected to have the form 

CW + (0,0,t) = Z x e- Emt + Z 2 e -JWCr-t) + Zg e -^ B * e -B B ,(T-t) 

+ Z 4 e - EBB{T - t] + Z 5 e^s'e'* + . . . , (11) 

where Egg is the energy of the two-baryon state, Eg is the ground-state baryon energy of 
positive parity, Eg> is the energy of the ground-state negative-parity baryon or meson-baryon 
scattering state and the ellipsis denotes contributions from higher excited states of the same 



quantum numbers. The exponent of the first term in eq. (11) is the primary object of our 
study and the other contributions are viewed as "pollutants". 

Since the interpolating operators for the baryons that give rise to the baryon-baryon 
correlation functions are individually projected with the positive-energy projectors, Y± = 
(1 ±7o)/2, the Z-factors of the forward and backward propagating states of the same energy 



in eq. (11) are not related, for example it is possible that Z$ ^ Z\. Indeed, since the sink 
separately projects the parity of each baryon (along with its momentum), the amplitude of 
the two backwards propagating baryons is suppressed by the lattice volume relative to the 
two forwards propagating baryons. This is a significant suppression, and we expect to see 
more energetic states, such as B'B', dominating the correlation function at large times near 
the boundary. 



C. Computational overview 

In the present study, we have employed a single ensemble of 20 3 x 128 gauge-field configu- 
rations generated with a rif = 2 + 1 flavor, anisotropic-clover quark action, with a spatial 
lattice spacing of b = 0.1227(8) fm and a pion mass m n ~ 390 MeV, that have been produced 
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by the Hadron Spectrum Collaboration [2U [25]. The technical details of the propagators 
computed on this ensemble are presented in Ref. [22] and we do not repeat them here. In the 
current calculation, the number of measurements has been increased to an average of 364 
randomly-distributed measurements on each of 1195 configurations (a total of ~ 0.435 x 10 6 
measurements). For correlators from each source point, two types of sink interpolating oper- 
ators are used [22] and the resulting correlation functions are referred to as "smeared-point" 
(SP) and "smeared-smeared" (SS). The measurements on each configuration are averaged, 
and then these averaged measurements are typically blocked (averaged) in sets of ten neigh- 
boring configurations (100 trajectories) to account for residual correlations (see Ref. [22] for 
a detailed study of correlations between different sources and configurations). 

The methods used to determine quantities of interest from the measured correlation func- 
tions are discussed in detail in Refs. [221 123], and we present no more than an overview here. 
One method of extraction is a direct analysis of "effective plots" . Linear combinations of the 
SP and SS correlation functions associated with each baryon and baryon-baryon state are 
formed to eliminate the contribution from excited states, as discussed in Ref. [22] . Effective 
mass plots are formed from these correlation functions and ratios of correlation functions 
to extract the energy and energy splitting, respectively 2 . Time intervals over which the 
effective mass appears constant are identified, and the energy is extracted from a correlated 
% 2 -squared minimization with the covariance matrix determined with the Jackknife or Boot- 
strap procedure. A fitting systematic uncertainty is assigned from the range of extracted 
energies determined by varying the location and size of the fitting interval over a reasonable 
range. A given energy-splitting and its uncertainty is converted into a value of g 2 and its 
uncertainty using eq. or a direct fit to an effective-g 2 plot is performed, which is then 
translated into q n cot 5(q n ) and its associated uncertainty using eq. ([8]). There are a num- 
ber of ways to perform this last stage of the analysis. One way is to use the Jackknife or 
Bootstrap procedure to determine the uncertainty in q n cot5(q n ) directly. However, this is 
complicated by the fact that S(x) in eq. ^ is a singular function. An alternate method 
is to determine the value of g 2 and its associated uncertainty from the two-baryon energy 
splitting from eq. ([6]), and then to propagate the central value and uncertainties through 
eq. ^ to determine q n cot S(q n ). We present results using the latter method. 

Results from this methodology are consistent with those from multi-exponential fits to the 
correlation functions and with various matrix-Prony [23] based analyses. In the following, 
we present results from a single analysis, ensuring that the systematic uncertainties are 
sufficient to maintain agreement with analyses using these other methods. 



2 Linear combination of SP and SS correlation functions are formed that extend the plateau of the ground 
state found from a matrix-Prony analysis [23 to earlier times. The effect of choosing a slightly different 
linear combination is taken into account by the systematic error. While we find that the best candidate 
for the ground state level is not ambiguous in any of the two-baryon correlation functions, there is always 
the possibility that the "true" ground state becomes dominant in the region where it is likely that there 
are significant contributions from backward-propagating states which contaminate the signal. Eliminating 
this possibility requires better statistical precision and/or measurements on lattices with a longer time 
extent. 
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III. SINGLE BARYONS 



The interaction between baryons is extracted from the difference between the energy-levels of 
the two baryon system in the lattice volume and the individual baryon masses. The masses 
of the baryons were extracted in a previous work [22], but it is useful to show the masses 
here, particularly due to the substantial increase in the number of measurements that have 
been performed. The baryon masses are extracted from correlated x 2 minimizing fits to the 
generalized effective mass plots (GEMP's) obtained from each correlation function, Ci(t), 
defined to be 



I 1 ( C ® 

tj ° g \C(t + tj) 



(12) 



The GEMP's from the single-baryon correlation functions are presented in fig. [TJ The 





FIG. 1: The single-baryon GEMP's (with tj = 3) resulting from the linear combination of SS 
and SP correlation functions that eliminates contributions from excited states. The fit values of 
the masses along with the statistical uncertainty, and the systematic and statistical uncertainties 
combined in quadrature, are shown. 



extracted baryon masses that are fit to the plateau regions of the GEMP's are shown along 
with the statistical uncertainty, and the systematic and statistical uncertainties combined in 
quadrature. The results of the fitting, along with the fitting intervals are given in Table |T] 

The statistical and fitting systematic uncertainties are below 0.2% for each of the single 
baryon states (although the uncertainty in the temporal lattice spacing leads to a larger 
uncertainty on the masses in physical units). The baryon signal-to-noise ratio is essentially 
independent of time in the fitting windows shown in Table [I] due to the nature of the sources 
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TABLE I: Extracted single hadron masses. A lattice spacing of b s = 0.1227 ± 0.0008 fm and an 
anisotropy factor of £ s = 3.500 ± 0.032 is used to convert from temporal lattice units (t.l.u.) to 
MeV. The first two uncertainties are the statistical and systematic uncertainty of the extraction 
in temporal lattice units, while the third uncertainty quoted for quantities in physical units is the 
combined lattice spacing and anisotropy uncertainty. 



Hadron 


M (t.l.u.) 


M (MeV) 


x/dof 


fitting 


; interval 


7T 


0.06936(12)(05) 


390.39(0.67)(0.28)(4.38) 


0.73 


21 


-» 41 


K 


0.097016(99)(33) 


546.06(0.56)(0.19)(6.13) 


1.01 


29 


-» 49 


N 


0.20682(34) (30) 


1164.1(1.9)(1.7)(13.1) 


1.16 


21 


-» 40 


A 


0.22246(27) (27) 


1252.1(1.5)(1.5)(14.1) 


0.97 


13 


38 


E 


0.22752(32) (29) 


1280.6(1. 7)(1.6)(14.3) 


1.46 


21 


-» 40 




0.24101 (27) (27) 


1356.5(1.5)(1.5)(15.2) 


1.06 


16 


-» 40 



that generate the correlation function, as discussed in Ref. [22]. It is only in these time 
intervals (for this value of tj) that reliable energy splittings between the two-baryon and the 
single-baryon masses can be constructed from the ratio of the correlation functions given in 
eq. @. 

IV. BARYON-BARYON INTERACTIONS 

The energy eigenstates in the finite lattice volume are classified by their global quantum 
numbers: baryon number, isospin, third component of isospin, strangeness, total momen- 
tum, and behavior under hyper-cubic transformations. Six quark operators that are simple 
products of three-quark baryon operators are used as sources for the baryon-baryon cor- 
relation functions. As a consequence, the baryon content of the interpolating operator is 
used to define the operator, e.g. nA( 3 Si), but this operator will, in principle, couple to all 
states in the volume with the quantum numbers B — 2, I — ^, I z = — |, s = — 1, and 
2s+1 Lj = 3 S\ + where the ellipses denote states with higher total angular momentum 
that also project onto the A\ irreducible representation of the cubic-group 3 . Both SS and 
SP correlation functions have been calculated for the nine baryon-baryon channels shown in 
Table OH 

If the calculations were performed on gauge-field configurations of infinite extent in the 
time-direction, so that only forward propagation could occur, some of the channels in Table [TT| 
could be analyzed by considering contributions from a single scattering channel, e.g. NN, 
nYr , as we expect a single, well-separated ground state for these quantum 
numbers. However, other channels may require a multi-channel analysis, e.g. nA, AA. The 
nA source will produce low-lying states in the lattice volume that are predominately linear 



The spatial dimensions of the gauge-field configurations that are used in this work are identical (i.e. 
isotropic), and as such the eigenstates of the QCD Hamiltonian can be classified with respect to their 
transformation properties under cubic transformations, H(3), a subgroup of the group of continuous three- 
dimensional rotations, 0(3). The two-baryon states that are calculated in this work all belong to the 
representation of H(3), corresponding to combinations of states with angular momentum L = 0, 4, 6, . . . . 



9 



TABLE II: Baryon-baryon channels examined in this work. 



1 ,h ATI Tl^l 


r 

_/ 


1 7~ 1 

\ 1 z\ 


Q 


fin \ 

PP ( ^o) 


1 


1 


A 

u 


np { bi) 


A 

u 


A 

u 


A 

u 


nA ('So) 


1 

2 


1 
2 


— 1 


nA ( 3 5i) 


1 
2 


1 
2 


-1 


nS- (%) 


3 
2 


3 
2 


-1 


nS- ( 3 Si) 


3 
2 


3 
2 


-1 


E-E- (%) 


2 


2 


-2 


AA (%) 








-2 


S-S- (!fib) 


1 


1 


-4 



combinations of the nA, nE° and pE~ two-baryon states. The AA source will produce low- 
lying states in the lattice volume that are predominately linear combinations of the AA, 
E ± '°E :F ' C) , and NE two-baryon states. 

A. Nucleon-nucleon interactions 

Perhaps the most studied and best understood of the two-hadron systems are the proton- 
proton and proton-neutron. At low energies, only two combinations of spin and isospin 
are possible, a spin-triplet isosinglet np ( 3 Si) and a spin-singlet isotriplet pp (}Sq). At the 
physical pion mass, the scattering lengths in these channels are unnaturally large and the 3 SI 
channel contains a shallow bound state, the deuteron, with a binding energy of ~ 2.22 MeV. 
These large scattering lengths and the shallow bound state arise because the coefficient of 
the momentum-independent four-nucleon operator in the low-energy effective field theory 
has a non-trivial ultraviolet fixed-point for the physical light-quark masses. An interesting 
line of investigation is the study of the scattering lengths as a function of the quark masses 
to ascertain the sensitivity of this fine-tuning to the QCD parameters (301 EH E2]- The fine 
tuning is not expected to persist away from the physical masses and we expect our present 
(unphysical) calculations to yield scattering lengths that are natural-sized. 

The GEMP obtained from the proton-proton correlation functions is shown in fig. [2j and 
fig. [3] shows the analogous plots for the neutron-proton correlation function. After the initial 
plateau region, the GEMPs show a slight downward fluctuation at time-slice ~ 29, which we 
believe is statistical in nature. 4 Figure [4] shows the effective |k| 2 plot 5 for both the proton- 
proton and neutron-proton channels. Both channels exhibit plateaus in |k| 2 . The plateau in 
the neutron-proton channel is consistent with zero, while the plateau in the proton-proton 
channel differs from zero at the ~ la-level. We conclude that at this value of the pion 



4 This fluctuation appears in both the proton-proton and neutron-proton GEMPs and in several of the 
other channels. The fact that the feature appears in multiple correlation functions is not a surprise as all 
the correlation functions are generated from the same light-quark and strange-quark propagators. 

5 For the presentation of the results of the calculation we use |k| 2 to denote q^, as the GEMPs do not 
isolate a particular energy-eigenvalue or eigenstate. 
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FIG. 2: The left panel is the proton-proton ( 1 5'o) GEMP with tj = 1, while the right panel shows 
the plateau region of the left panel. The band in the right panel and the upper line in the left 
panel correspond to 2Mjy, while the lower two lines in the left panel correspond to —2Mn and 
— 2(Afjv + m n ), respectively. 



mass, the interactions between nucleons produce a small scattering length in both channels 
compared to the naive estimate of m~ l ~ 0.5 fm. The results of the analysis are shown in 
Table III. Extracted lattice quantities are converted to physical units using b s = 0.1227(8) 
fm; here the lattice spacing uncertainty is sub-dominant. Motivated by the fact that the pion 



TABLE III: Results for the pp ( 1 5'o) and np ( 3 <Si) channels. 
Process |k| 2 /m^ AE (MeV) — 1/pcot 5 (fm) x 2 /dof fitting interval 



pp 0.030(13)(20) 3.9(1. 7)(2.6) 0.118+°;°^°;°^ 1.6 20 ^32 

np 0.012(20)(33) 1.6(2.6)(4.3) 0.052+°;°^°'JJ 1.96 19^32 



mass dictates the range of the interaction between nucleons, we have shown m n /p cot 5(p) 
as a function of |k| 2 /m 2 in fig. [HJ 

A summary of lattice QCD calculations of NN scattering is shown in fig. [6] Since the 
momenta at which the phase shifts are measured are small, we present these results as 
scattering lengths, implicitly assuming that higher order coefficients in the effective range 
expansion, eq. (10), are natural sized. The results calculated in this work are consistent with 
those that we obtained using mixed-action lattice QCD [19J. It is interesting to note that 
the results of quenched calculations [29] yield scattering lengths that are consistent within 
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FIG. 3: The left panel is the neutron-proton ( 3 Si) GEMP with tj = 1, while the right panel shows 
the plateau region of the left panel. The band in the right panel and the upper line in the left 
panel correspond to 2Mjy, while the lower two lines in the left panel correspond to —2Mn and 
— 2(Mjv + rn n ), respectively. 
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FIG. 4: The left panel is the effective |k| 2 plot for the proton-proton ( 1 5'o) channel with tj 
and the fit to the plateau. The right panel is for the neutron-proton ( 3 Si) channel. 



uncertainties with the fully-dynamical rij = 2 + 1 values. 
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FIG. 5: The inverse of the real part of the inverse scattering amplitude normalized to the pion mass 
as a function of the squared-momentum in the center-of-mass normalized to the pion mass. The 
dotted curve corresponds to the inverse "S-function" , defined in eq. Q, from which (fccot^) -1 
is determined from |k| 2 . The inner uncertainty of each data point is statistical and the outer 
uncertainty is the statistical and systematic uncertainty combined in quadrature. 
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FIG. 6: A compilation of the scattering lengths for NN scattering in the 1 Sq (left panel) and 3 S*i 
(right panel) calculated in lattice QCD and quenched lattice QCD. The data are from Refs. [28 , 
|19j . |29] and the current work. The vertical dashed line corresponds to the physical pion mass. 
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B. Hyperon-nucleon interactions (s = — 1) 



Lattice QCD calculations of the YN interactions are of greater phenomenological impor- 
tance than those of NN interactions because of the limited experimental access to hyperon 
systems in the laboratory and the possible role of hyperons in nuclear astrophysics. We have 
calculated the energy-eigenvalues of systems with the quantum numbers of nA and nS~ in 
both spin channels. 



1. nS interactions (I = |J 

Interpolating operators with the quantum numbers of nS~ in either the 1 5'o or 3 Si chan- 
nels will couple to energy eigenstates in the lattice volume with strangeness s = — 1 and 
isospin of / = |, and the spectrum is not expected to have more than one low-lying ground 
state. We expect to be able to describe these systems with a single elastic-scattering channel 
for the lattice volumes we are working in as these are the only two-baryon states comprised 
of octet baryons that have these quantum numbers. 

The GEMP for the nYi~ ( 1 5'o) is shown in fig. [7] and exhibits a clear plateau, as does 
the GEMP for the ( 3 Si) channel that is shown in fig. [8j Figure [9] shows the effective 




FIG. 7: The left panel is the ( 1 5'o) GEMP with tj = 1, while the right panel shows the 

plateau region of the left panel. The band in the right panel and the upper line in the left panel 
correspond to My + Mn, while the lower two lines in the left panel correspond to —{My, + M^) 
and —{My, + Mn + 2m 7r ), respectively. 

|k| 2 plot for both the nX~ ( 1 5'o) and nS~ ( 3 Si) channels. The results of fitting the clear 



plateaus that are observed in both channels are shown in Table |IV[ and fig. 10 shows the 



results presented in Table IV normalized to the pion mass. The nE ( 1 5'o) channel is observed 



to have an interaction of natural size, and the energy that is measured in the calculation 
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FIG. 8: The left panel is the nS~ ( 3 Si) GEMP with tj = 1, while the right panel shows the 
plateau region of the left panel. The band in the right panel and the upper line in the left panel 
correspond to Ms + M^v, while the lower two lines in the left panel correspond to —{My, + M^r) 
and —(Ms + Mn + 277%), respectively. 




FIG. 9: The left panel is the effective |k| 2 plot for the nS ( 1 S'o) channel with t j = 1 and the fit 
to the plateau. The right panel is for the raS - ( 3 S\) channel. 

lies within the regime of applicability of the effective-range expansion. In contrast, the 
interaction in the n£~ ( 3 5'i) channel is seen to be large, \m n /p cot 5(p)\ ~ 3, and is well 
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TABLE IV: Results for the strangeness s = 



— 1 hyperon-nucleon channels. 



Process |k| 2 /m 2 AE (MeV) — 1/pcot 5 (fm) x 2 /dof fitting interval 



2.06 20 -> 35 

0.91 19 -> 36 



1.69 20 -> 35 

1.05 20 -> 35 




FIG. 10: The inverse of the real part of the inverse scattering amplitude normalized to the pion 
mass as a function of the squared-momentum in the center-of-mass of the baryons normalized to 
the pion mass. The dotted curve corresponds to the inverse "S-function" , defined in eq. (pi), from 
which (kcot5)~ 1 is determined from |k| 2 . The inner uncertainty of each data point is statistical 
and the outer uncertainty is the statistical and systematic uncertainty combined in quadrature. 

outside the regime of applicability of the effective range expansion. At this momenta, the 
n£~ ( 3 5i) is strongly interacting in a way that is consistent with an attractive interaction 
that supports a bound-state (which we find no direct evidence for in this calculation), or 
a repulsive interaction of unnaturally large range. These two scenarios cannot be resolved 



nE- ('So) 0.122(12)(19) 15.3(1.5)(2.3) 0.361±™|™ 
nE~ ( 3 Sr) 0.551(17)(19) 67.9(2.1)(2.3) 1.47i°;^°^ 



nk (%) 0.093(12)(19) 11.8(1.6)(2.3) 0.297^™ 
nA ( 3 5x) 0.094(15)(15) 11.9(1.9)(1.9) ^mt^-lZ 
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with calculations in a single volume, but ongoing calculations in different volumes will resolve 
this ambiguity. This result is perhaps the most important result of this present work. One 
concludes from this calculation that the nY~ interaction is strongly spin-dependent. 



2. The coupled NA- NE channel (I = \) 

The strangeness s = — 1, isospin-^ energy-eigenstates in the lattice volume will, in gen- 
eral, couple to both the NA and NY channels. Therefore, the calculated NA correlation 
functions will receive contributions from all such eigenstates in the volume. However, as 
M s — Ma = 0.0051 t.l.u. is a significant mass splitting, we expect that a single-channel 
analysis is applicable in this system and proceed accordingly. Nevertheless, calculation of 
the isospin-^ NY correlation function and the crossed correlation function resulting from a 
iVA source and a NY sink would likely improve this analysis. 

The GEMPs for the nA ( 1 S'o) and nA ( 3 Si) correlation functions are shown in fig. 



and 12, respectively, and in both cases a clear plateau is visible. Figure 13 shows the 
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FIG. 11: The left panel is the nA ( 1 5'o) GEMP with tj = 1, while the right panel shows the 
plateau region of the left panel. The band in the right panel and the upper line in the left panel 
correspond to Ma + Mjy, while the lower two lines in the left panel correspond to —(Ma + Mat) 
and —(Ma + Mjy + 2??^), respectively. 

effective |k| 2 plot for both the nA ( 1 5'o) and nA ( 3 Si) channels and clear plateaus are again 
observed in both channels. The results of fitting to the plateau region of the effective |k| 2 



plots are reported in Table |IV[ and are displayed in fig. [10} The energy-splittings are both 
found to be AE ~ 0.002 t.l.u, and are therefore smaller than the expected splitting between 
the energy eigenstates resulting from the NA-NY mixing (M^ — Ma = 0.0051 t.l.u.). It 
therefore seems likely that, a posteriori, the single-channel analysis used here is applicable. 
The extracted squared-momenta and hence scattering amplitudes in the two spin channels 
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FIG. 12: The left panel is the nA ( 3 Si) GEMP with tj = 1, while the right panel shows the 
plateau region of the left panel. The band in the right panel and the upper line in the left panel 
correspond to Ma + Ma, while the lower two lines in the left panel correspond to —(Ma + Ma) 
and —(Ma + Mat + 2m n ), respectively. 
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FIG. 13: The left panel is the effective |k| 2 plot for the nA ( 1 5'o) channel with t j = 1 and the fit 
to the plateau. The right panel is for the nA ( 3 Si) channel. 

are the same within uncertainties. This indicates that the spin-dependent interactions in 
these channels are very small. If the extracted states are predominately iVA, then this result 
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TABLE V: A compilation of results for strangeness = — 1 hyperon-nucleon scattering from lattice 
QCD. The columns labeled as "Valence" and "Sea" list the action used for the valence and sea 
quarks. A dash indicates a quenched calculation. 



Process 


m n (MeV) 


|k| (MeV) 


— 1/pcot 5 (fm) 


Valence 


Sea 


Reference 


nA 


So 


296(3) 


50(26) i 


-0.11(7) 


Clover 


Clover 


[33] 


nA 


'S 


354(6) 


255(26) 


1.04(28) 


DW 


Staggered 


[20 


nA 


1 So 


390.4(4.4) 


119(14) 


0.297(76) 


Clover 


Clover 


present work 


nA 


'So 


465(1) 


22(3) i 


-0.09(3) 


Clover 




|3"51 

1 J 


nA 


'So 


493(8) 


197(24) 


0.63(12) 


DW 


Staggered 


120 


nA 


'So 


514(1) 


17(3) i 


-0.07(3) 


Clover 




133 

L 


nA 


3 Si 


296(3) 


40(24) i 


-0.07(7) 


Clover 


Clover 


1 33 1 

L J 


nA 


3 s 1 


354(6) 


168(65) 


0.50(27) 


DW 


Staggered 


1201 


nA 


3 s 1 


390.4(4.4) 


119(14) 


0.299(49) 


Clover 


Clover 


present work 


nA 


3 S 1 


465(1) 


24(3) i 


-0.11(3) 


Clover 




[33] 


nA 


3 s 1 


514(1) 


20(2) i 


-0.09(2) 


Clover 




m 


nE~ 


'So 


390.4(4.4) 


136(12) 


0.361(46) 


Clover 


Clover 


present work 


n£~ 


'So 


493(8) 


179(30) 


0.57(14) 


DW 


Staggered 


[20] 


n£~ 


3 s 1 


390.4(4.4) 


289.8(6.8) 


1.47(16) 


Clover 


Clover 


present work 


n£~ 


3 Si 


493(8) 


261(37) 


1.19(53) 


DW 


Staggered 


m 


nS~ 


3 Si 


592(10) 


226(30) 


0.85(22) 


DW 


Staggered 


[20. 



is expected because of the fact that the long-range spin-dependent interaction resulting from 
OPE is absent (as the A is an isosinglet). Further, the extracted squared-momenta are 
consistent with that found in the nE - ( 1 5'o) channel. 

3. Compilation of measurements 

After the pioneering quenched calculations of YN scattering by Fukugita et al [28j, and 
the first fully-dynamical rif = 2 + 1 calculations [20], more refined quenched calculations 
have been performed along with one rif = 2 + 1 calculation [33]. All of the results for s = — 1 
YN scattering that have been obtained from lattice QCD calculations are shown in Table [Vj 
As the measurements have been performed at different pion masses and in different lattice 
volumes, resulting in different center-of-mass energies, it is difficult to present these results 
in a single diagram. The present measurements and the mixed-action measurements [20] 
were both on lattices with spatial volumes of V ~ (2.5 fm) 3 , and with a lattice spacing of 
b ~ 0.125 fm and so a direct comparison is possible for this subset of measurements. 

C. Hyperon-hyperon interactions (s = —2) 

Lattice QCD calculations of hyperon-hyperon interactions are important as they can provide 
guidance to the experimental programs in hyper-nuclear physics. They can also improve 
upon the current understanding of the stability of the core of supernovae if it becomes 
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energetically favorable to have strange baryons present. To this end we have calculated the 
correlation functions resulting from AA and E~E~ sources and sinks. The E~E~ channel 
has I z = —2 and, as such, a single, isolated ground-state is expected which can be used to 
determine the E~E~ scattering phase shift with a single-channel analysis. In contrast, the 
correlation function produced by the AA source is expected to exhibit two nearly degenerate 
states as 2M A = 0.44492 t.l.u. and M s + M N = 0.44783 t.l.u.. Further, the E+E" state 
is likely to be close by, 2M% = 0.45504 t.l.u.. An operator of the form E~p would help 
to resolve the states, but we have not calculated this correlation function, nor the possible 
mixed correlation functions. Further, we have not explored the isospin-1 channels, such 
as the NE — AS coupled channels, but we would expect to find very-closely spaced energy 
eigenstates for the pion mass used in this calculation, making clean extractions very difficult. 
In particular, M E -+M n = 0.449 t.l.u while M A + M S - = 0.452 t.l.u. Quenched calculations 
of the pE° scattering length have been presented in Ref . [34J . 



1. A A interactions (1 = 0) 



The low-lying eigenstates that couple to a AA source and sink are, in principle, linear 
combinations of all two-baryon states with s = —2 and 1 = 0, namely, AA, EE, NE and 
their excitations. The GEMP for the AA channel (for a particular combination of SP and SS 




FIG. 14: The left panel is the AA GEMP with tj = 3, while the right panel shows the plateau 
region of the left panel. The band in the right panel and the upper line in the left panel correspond 
to 2Ma, while the lower two lines in the left panel correspond to —2M\ and — 2(Ma + m n ), 
respectively. 



correlation functions) is shown in fig. 14 , in which a clear plateau is observed. The effective 
|k| 2 plot shown in fig. 15 also shows a clear plateau that is negative shifted in energy, 



unambiguously indicating an attractive interaction. The energy of this state is the lowest, 
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FIG. 15: The effective |k| 2 plot for the AA channel with t j = 3 and the fit to the plateau. 



and the plateau is the longest, for any combination of SS and SP correlation functions that 
could be constructed. Other combinations of correlation functions suggest states at higher 
energies, consistent with expectations of nearby states, however, with the current data, no 
definitive statements can be made. The splittings between the asymptotic states, AA and 
iVS are such that at measured (very small) AA center-of-mass momentum a single-channel 
analysis can be performed to define the AA elastic scattering amplitude, the results of which 



are shown in Table VI, and shown graphically in fig. 16 



TABLE VI: Results for the strangeness = — 2 and strangeness = —4 hyperon-hyperon channels. 



Process 


\k\ 2 /ml 


AE (MeV) 


— 1/pcot 5 (fm) 


X 2 /dof 


fitting 


interval 


AA 


-0.033(09)(11) 


-4.1(1.2)(1.4) 


n 1OQ+0.062+0-072 
u - ioo -0.072_o.085 


1.38 


20 


-> 35 


£-£- 


0.215(11)(18) 


25.5(1.3)(2.1) 


r:q 4 +0. 019+0-032 
U,OM -0. 019-0. 032 


2.6 


20 


-> 36 




0.0247(94)(77) 


2.8(1. 1)(0.9) 


n 1 m +0.032+0-026 
u - lul -0. 036_o. 029 


1.56 


21 


-> 34 



21 




0.5 



\k\ 2 /ml 



FIG. 16: The inverse of the real part of the inverse scattering amplitude normalized to the pion 
mass as a function of the squared-momentum in the center-of-mass of the baryons normalized to 
the pion mass. The dotted curve corresponds to the inverse "S-function" , defined in eq. from 
which (/ccot^) -1 is determined from |k| 2 . The inner uncertainty of each data point is statistical 
and the outer uncertainty is the statistical and systematic uncertainty combined in quadrature. 



The AA channel is the only channel in which a negatively-shifted energy splitting is 
observed. However, without performing measurements on additional lattice volumes, it is 
presently not possible to determine if this negative energy shift indicates the presence of 
a bound state (the infamous H-dibaryon [35J) or if it is simply a continuum state that is 
negatively shifted due to an attractive interaction. The location of the state on the "S- 
function" curve suggests that it is fact a continuum state [8] , but further measurements are 
required to properly explore this exciting possibility. 



2. £ £ interactions (I = 2) 



Due to its quantum numbers, s = —2 and 1 = 2, the channel is not expected to 

have other states near the ground state in the lattice volume as there are no other states 
comprised of two octet baryons with these quantum numbers. The GEMP associated with 
the S~S~ correlation function is shown in fig. 17 and the resulting effective |k| 2 plot is shown 
18 The results obtained by fitting to the effective |k| 2 plot are shown in Table VI 



in fig. 
above. 



and presented graphically in fig. 16 



due to the downward fluctuation near time-slice t 



The somewhat large value of x 2 /dof 
29 in fig 



2.6 is 

which we assume to be 



statistical in nature. A smaller fitting interval would yield approximately the same scattering 
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FIG. 17: The left panel is the GEMP with tj = 1, while the right panel shows the plateau 

region of the left panel. The band in the right panel and the upper line in the left panel correspond 
to 2Ms, while the lower two lines in the left panel correspond to — 2M^ and — 2(Afs + m n ), 
respectively. 

phase shift, but with a substantially smaller x 2 /dof. 



D. Hyperon-hyperon interactions (s = —4) 

Baryon-baryon interactions in the strangeness -4 sector have no obvious phenomenological 
implications. However, it is found that systems containing heavier quarks, such as the Q or 
S, have better-behaved correlation functions in lattice QCD calculations. As lattice QCD 
calculations of hadronic interactions are still in their infancy, it is useful to explore such 
systems to better understand aspects of the methodology that we employ. 

The only s = —4 systems of two octet-baryons are the I = 0, 1 HS combinations and 
here we focus on / = 1 (H - E~). The low- lying states in the lattice volume that couple to 
the S~S~ interpolating operator are expected to be describable in terms of a single-channel 
elastic-scattering amplitude as there are no other states composed of two octet baryons that 



can couple to them. The GEMPs for this channel are shown in fig. 19, and the resulting 



effective |k| plot is shown in fig. 20 The results obtained by fitting to the effective |k| plot 



are shown in Table |VI[ and presented graphically in fig. [16j normalized to the pion mass. 
The effective |k| 2 plot shows the downward fluctuation seen in other correlators at time-slice 
t = 29, and the resulting fit is consistent with zero at the 2<r-level. We conclude that the 
interactions at this value of the pion mass are quite weak. 
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FIG. 18: The effective |k| 2 plot for the channel with tj = 1 and the fit to the plateau. 

V. STATISTICAL SCALING AND NOISE IN CORRELATION FUNCTIONS 

The precision of lattice QCD calculations of any quantity is limited by the statistical noise 
in the relevant correlation functions. Until recently [22} [23] , the lore has been (based on gen- 
eral arguments by Lepage [18]) that the signal-to-noise ratios in (multi-)baryon correlation 
functions degrade exponentially with time and with the number of baryons in the system. In 
Refs. [22l [23], we showed that, while this behavior is observed at large times, at intermediate 
times the signal-to-noise ratio does not degrade exponentially, and in fact it is found to be 
independent of time for a significant number of time slices (the "Golden Window") for the 
sources that are used. Further, in this window of time-slices, the signal-to-noise ratio is 
essentially independent of the number of baryons in the system. 

After reviewing Lepage's arguments regarding the general behavior of signal-to-noise 
ratio, and their generalization to the temporal boundary conditions that are used in the 
present work, a thorough exploration of the noise in the various correlators of the two- 
baryon sector is presented. 

A. Aspects of the noise correlation functions and the signal-to-noise ratio 

On gauge-field configurations with infinite temporal extent, correlation functions of one 
or more baryons exhibit statistical noise at large-times that increases exponentially with 
Euclidean time [18]. In the case of a source that has the quantum numbers of a single 
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FIG. 19: The left panel is the GEMP with tj = 1, while the right panel shows the plateau 

region of the left panel. The band in the right panel and the upper line in the left panel correspond 
to 2Ms, while the lower two lines in the left panel correspond to — 2M= and — 2(M= + m n ), 
respectively. 



positive parity nucleoli, the correlation function has the form 

r^/n n\ lr>\ '7 „-M N t 



(t)) = J2 T + <°l N a (x,t)W(0,0) |0) - Z N e 



(13) 



where iV a (x, t) is an interpolating field (composed of three quark operators) that has non- 
vanishing overlap with the nucleon, r + is a positive energy projector, and the angle brackets 
indicate statistical averaging over measurements on an ensemble of configurations. The 
variance of this correlation function is given by 



Na 2 ~ (0]v(t)M*)> - (W)) 2 

= J2 T + Tl f <°l N a (x,t)N\y,t)N~<(0,0)N\o,0) |0> 



y _ a -2M N t 7 2 -2M N t , r? „-3nM 



t— >oo 



v.v ~~ ^.x' + e am "* 1 + ... — > Zz-n e 3mnt , (14) 



where all interaction energies have been neglected, and N is the number of (independent) 
measurements (distinct from the nucleon field operator N). At large times, the noise-to- 
signal ratio consequently behaves as [18] 



a 

x 



a(t) 



1 



More generally, for a system of A nucleons, the noise-to-signal ratio behaves as 
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FIG. 20: The effective |k| 2 plot for the 



channel with t j = 3 and the fit to the plateau. 



at large times. 

As we discussed in Ref. [23], the various "Z-factors", such as Z 37T , depend upon the 
details of the sources and sinks interpolators that are used. For the present calculations, the 
projection onto zero-momentum final state nucleons, introduces a 1 / -\/Volume suppression 
of the amplitudes of the various components (except for NN) in addition to color and 
spin rearrangement suppressions that exists independent of the spatial structure of the 
source. As a consequence, an interval of time slices exists at short times (the "Golden 
Window") in which the variance of the correlation function is dominated by the terms in 



eq. (14) that behave as ~ e~ N . In this window, the signal-to-noise ratio of the single 
baryon correlation function is independent of time. Further, the signal-to-noise ratio does 
not degrade exponentially faster in multi-baryon correlation functions than in single-baryon 
correlation functions in the "Golden Window" [23] . 

The finite temporal extent introduces backward propagating states (thermal states) into 
the correlation functions which lead to exponentially worse signal-to-noise ratios at large 
times [221 123] • These contributions are suppressed by at least exp(m 7r T), however, in the 
present work (where m^T ~ 9), these effects cause complications. We note that the impact 
of these states can be mitigated by working at larger temporal extents and exponentially 
large computational resources are not required to remove this effect. 

With the high statistics that have been accumulated in the present work, the behavior 
of the signal-to-noise ratio can be carefully examined. It is useful to form the effective 
noise-to-signal plot [22], in analogy with the GEMPs. On each time slice, the quantity 



S(t) 



x(t) 



(17) 



26 



is formed, from which the energy governing the exponential behavior (the signal-to-noise 
energy-scale) can be extracted via 



"i^Ber) • (18) 

For a correlation function that is dominated by a single state with a corresponding vari- 
ance correlation function dominated by a single energy scale, the quantity E s (t;tj) will be 
independent of both t and tj. 



B. Measured signal-to- noise ratios in the one and two nucleon sectors 



In the single nucleon sector, we expect that the energy scales E s ~ 0, M 



jv — §777.71-, and others, 

contribute to the signal-to-noise ratio. At times when the nucleon correlation function is 
in the ground state, and the variance correlation function is dominated by the nucleon- 
antinucleon state, E s = should dominate the signal-to-noise ratio. At large times the 
variance correlation function is dominated by the 3-pion state and E s = — ^m^ should 
dominate. This is modified by the finite temporal direction [22] as the hadrons produced by 
the sources of the correlation function and the variance correlation function can propagate 
forward and backward in time. The measured energy scale of the signal-to-noise ratio of the 
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FIG. 21: The energy scale of the signal-to-noise ratio in the nucleon (left panel) and proton-proton 



(right panel) correlation functions, as defined in eq. (18), with tj = 6. The horizontal lines in 



the left panel correspond to E s = 0, Mn — and — ^m n , while those in the right panel 
correspond to E s = 0, Mjy — \m % , Mjy — \m 1T , 2Mn — Sm^, 2Mn — m^. 



single nucleon correlation function is shown in fig. [2TJ It exhibits behavior that is consistent 
with expectations, and exceeds the long-time behavior expected from the Lepage argument 
at approximately time-slice t = 50 due to the temporal boundary conditions. 
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On configurations with infinite temporal extent, the proton-proton correlation function 
is of the form (neglecting interactions between the hadrons) 



(W*)> = E V + PP (°l ^CM^Cy, t)iV>,0)iV P (O,0) |0) 

- e- 2M ^ + . . . , (19) 
and the variance correlation function has the form 



Na 2 ~ (e]f N {t)e NN (t)) - (9 NN (t)) 2 

= J2 r^r? 7Ct (0| N a (x,t)N p (y,t)N\z,t)N v (w,t) x 

x,y,z,w 

^(0,0)^(0,0)^(0,0)^(0,0) |0) - (9 NN (t)) 2 

, v p -iM N t 7 2 -AMjxt , 7 _ -(2M JV +3m^)t , 7 „-6m,t , 

-> Z 67r e -6 " 1 **. (20) 

Therefore, we anticipate finding energy scales of approximately i? s = 0, — \m. K and 
2Mat — 3m7r in the signal-to-noise ratio on gauge-field configurations of infinite temporal 
extent. The temporal boundary conditions imposed in the present calculation introduce 
additional energy scales due to the backward propagating states. 



Figure 22 shows the energy scale associated with the signal-to-noise ratio for the ratio of 
correlation functions that provides the energy splitting between two interacting protons and 
two isolated protons from which the pcot<5(p) is extracted. It is clear that the energy scale 
of the energy splitting is significantly less than for the individual energies, and is consistent 
with zero throughout much of the Golden Window of time slices. This indicates that the 
signal-to-noise ratio associated with the energy splitting in the proton-proton sector, and 
hence the scattering parameters and bound-state energies, are time independent for the 
sources and sinks used here, and therefore do not degrade exponentially with time. This is 
an exceptionally important result, as it means that the extraction of NN, and more generally, 
multi-nucleon interactions, does not require an exponentially large number of measurements 
for each relevant correlation function. Further exploration and discussion of this point can 
be found in our work on three baryon systems [23] . 



C. Measured signal-to- noise ratios in the NT, sector 

The discussion of the behavior of the signal-to-noise ratio in the YN sector parallels that 



in the multi-nucleon sector in section VB An important difference is the presence of the 



strange quark, and the associated strange hadrons. The lowest energy scale contributing to 
the signal-to-noise ratio (beyond E s = 0) is E s = — ^rn n , and so it is expected that the 
degradation of the signal-to-noise ratio will be similar to that found in the NN correlations 
functions. It is found that the associated energy scale in the nE - ( 3 *Si) channel, as shown in 



fig. 23l is somewhat less than Mn — Im^ over a number of time-slices, and starts to exceed 



this value for time-slices greater than t> 20. 

The energy scale associated with the signal-to-noise ratio in the difference in energy 



between iVE interacting in the Si-channel and Mm + My, is shown in fig. 24. The energy 



scale seems to be somewhat larger than in the NN sector, and is non zero throughout the 
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FIG. 22: The energy scale of the signal-to-noise ratio, as defined in eq. (18), in the ratio of 



correlation functions that produces the shift in energy between two interacting protons and two 
isolati 
3m n . 



isolated protons, with tj = 6. The horizontal lines correspond to E s = 0, Mjy — \m lx and 2Mjy — 



Golden Window. The interaction is strong in this channel, and therefore a non-zero value 
of E s in the plateau region is not surprising. 



D. Measured signal-to-noise ratios in the one and two 3 sectors 

The signal-to-noise ratio in the EE sector is noticeably better than in the NN and the YN 
sectors. The energy scale associated with the H and E~E~ correlation functions are shown 



in fig. [25} The lowest energy scale contributing to the signal-to-noise ratio in the single H 
correlation function (beyond E s = 0) is E s = M= — mx — \m v . It is clear from the left panel 
in fig. 25 that the scale is much lower than this until time-slice t ~ 50. The situation is 
similar in the E~E~ correlation function in which the energy scale associated with the signal- 
to-noise ratio is found to be much less than the anticipated E s = M= — — \^r l (beyond 
E s = 0) until time-slice t ~ 40. Further, as shown in fig. 26, the energy scale associated 
with the energy splitting between the ETET state and two isolated S~'s is very small, and 
consistent with zero, for many time-slices below t< 35, and increases slowly beyond this. 

It is likely that the improved signal-to-noise behavior in the EE sector is due to a reduced 
overlap of the source onto the multi-meson intermediate states in the variance correlation 
function compared to purely baryonic intermediate states. Such a reduction is expected 
based on the fact that the volume occupied by multiple S's is smaller than that of multiple 
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FIG. 23: The energy scale of the signal-to-noise ratio in the £ (left panel) and the nS~ ( 3 Si) (right 
panel) correlation functions, as defined in eq. (18), with tj = 6. The horizontal lines in the left 
panel correspond to E s = 0, Ms — m# — \m^^ Ms — mx + and Ms — m^, while those in the 

3™ Jlyf ™ 1 



right panel correspond to E s 
M N + M s - m K - 2m n . 



0, M N - ±m w , M E - m K 



\M N + ±M S - \m K 



m 1 



and 



nucleons, and serves to extend the Golden Window beyond its range in nucleon correlation 
functions. 



E. Scaling of correlation functions, energy levels and scattering phase shifts 



In generating the 0.435 x 10 6 measurements on this ensemble of gauge-field configurations, 
we have performed an average of 364 measurements on each of the 1195 configurations. The 
scaling of the statistical and (fitting) systematic uncertainties in the single-hadron masses as 
a function of the number of measurements and number of gauge-field configurations in this 
ensemble was detailed in Ref. [22] • While the pion mass extraction was found to saturate as 
the number of measurements per configuration increased, the single baryon mass extractions 
did not saturate and scaled in a way that is approximately consistent with each measurement 
on the configuration being statistically independent. It is important to determine the scaling 
of uncertainties associated with scattering parameters determined in this lattice QCD calcu- 
lation because this scaling dictates the distribution of computational resources between the 
production of gauge-field configurations and the measurements performed per configuration. 
The AA channel provides a clean illustration of the scaling that is observed in the two-baryon 
sector. The fit to the blue circles in fig. 27 shows the scaling of the statistical uncertainty 
of the extracted value of q\ from the lowest level in the AA correlation function. The points 
correspond to the inverse of the variance of function of the number of sources per 

configuration, A^ rc , on 1155 gauge field configurations. The straight line corresponds to the 
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FIG. 24: The energy scale of the signal-to-noise ratio, as defined in eq. (18), in the ratio of 



correlation functions that produces the shift in energy between interacting S's and neutrons in the 
3 Si-channel and isolated S's and neutrons, with tj = 6. The horizontal lines correspond to E s = 0, 
Mat — nm-n and Ms — rriK — hm^. 



fit 1/er 2 = A iV s " c , with a = 0.94 ± 0.04. The fit to the purple squares in fig. 27 shows the 
scaling of the statistical and fitting-systematic uncertainties of the extracted value of q$ from 
the lowest level in the AA correlation function combined in quadrature. The straight line fit 
gives a = 1.11 ± 0.02. It is clear that both the statistical, and the combined statistical and 
systematic uncertainties, are scaling as ~ 1/ y/N^, consistent with statistically independent 
measurements even up to approximately 400 measurements per configuration. 

To emphasize the impact of the uncertainties on the extracted scattering parameters, and 
in particular, the need for high-statistics measurements of the scattering processes, fig. 28 
shows the extracted values of (jocot^)" 1 for AA scattering versus the extracted value of 
|k| 2 for different numbers of measurements per configuration, each with 1155 gauge-field 
configurations. With just 10 measurements per configuration the uncertainty in pcot5 
is large enough so that it is not possible to determine if the interaction is attractive or 
repulsive. This remains the case even for 100 measurements per configuration. It is only 
when the number of measurements per configuration approaches ~ 400 that the interaction 
can be determined to be attractive, but only with ~ 2a significance. 

Within uncertainties, the same scaling behavior is observed in all of the two-baryon 
channels. It appears that further measurements could be performed on this ensemble of 
gauge-field configurations that would continue to reduce the uncertainties in the two-baryon 
correlation functions. This statement is also valid for the single-baryon correlation functions 
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FIG. 25: The energy scale of the signal-to-noise ratio in the 3 (left panel) and (right panel) 

6. The horizontal lines in the left panel 
rriK + and Ms + mx 



correlation functions, as defined in eq. (18), with tj 
correspond to E s = 0, Ms — mx ' 



km^, while those in 



the right panel correspond to E s = 0, Mg — — lm„, M= — m# + im„, Ms + mj- — and 



2 (M= 



|m^) 



which do not show signs of saturation. 



VI. DISCUSSION 



We have calculated nucleon-nucleon, hyperon-nucleon and hyperon-hyperon interactions, 
with a high-statistics lattice QCD calculation on anisotropic improved-clover gauge-field 
configurations at a pion mass of m n ~ 390 MeV. A summary of the scattering information 
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that has been extracted from the measurements is presented in fig. 29 and in fig 

The phase shifts that we have obtained in the NN sector are small, and essentially con- 
sistent with zero. As a result we have been able to set a tight limit on the scattering 
lengths in both the 3 Si and 1 S'o channels (without extrapolation in |k| 2 ), as shown in fig. [6j 
These limits are consistent with our previous calculations of the scattering lengths in these 
channels |19j . but significantly more precise. 

Precise measurements in the YN sector (strangeness = — 1) have been obtained. The 
interaction in the n£~ ( 3 5'i) channel is found to be strong, but we are unable to determine 
if the interaction is attractive or repulsive until further measurements in different lattice 
volumes are performed. Such calculations are in progress. The measured momentum is far 
outside of the region for which an effective range expansion is convergent (|k| < m n /2), 
and consequently arguments concerning the naturalness or unnaturalness of the scattering 
amplitude are not possible. In contrast, the measured momentum in the n£ - ( 1 5'o) channel 
is within the region for which an effective-range expansion is convergent and the scattering 
length and effective range are either of natural size, or there are strong cancellations in the 
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FIG. 26: The energy scale of the signal-to-noise ratio, as defined in eq. (18), in the ratio of 
correlation functions that produces the shift in energy between interacting 
with tj = 6. The horizontal lines correspond to E s = and Mg — mx 



's and isolated S's, 



2 m n- 



effective range expansion. It is clear that the nE~ interactions are strongly spin dependent, 
as is expected from the long-distance contribution from one-pion exchange. On the other 
hand, the 1 5'o and 3 Si nA energy shifts are very similar, and as such we conclude that the 
interactions are essentially spin independent, as would be expected from channels without 
one-pion exchange. 

The measurement of a negatively-shifted energy level in the AA channel (strangeness 
= —2) indicates (at the statistical precision of the measurement) that the AA interaction is 
attractive. This is the only baryon-baryon channel for which we have measured a negative 
energy shift. This is an exciting measurement as it confirms that the channel in which 
the H-dibaryon [35] would arise is attractive. The present measurement suggest that the 
state does not correspond to a bound state in the infinite-volume limit at this pion mass 
(77^ ~ 390 MeV), but one can readily imagine that a bound state could arise at a lighter 
pion mass. 

The present work clearly demonstrates that, with sufficient computational resources, 
lattice QCD can be used to extract baryon-baryon scattering amplitudes as a function of 
momentum, and hence constrain the interactions between baryons. As has been discussed 
extensively in the literature [H |27j, it is not possible to directly extract the hadron-hadron 
potential (unless one or more of the quarks in each hadron is infinitely heavy), but effective 
interactions that reproduce the measured scattering amplitudes can be constructed and 
used in the calculation of other quantities of interest, in the same way that the modern NN 
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FIG. 27: The scaling of the statistical uncertainty (blue circles), and the statistical and fitting 
systematic uncertainties combined in quadrature (purple squares), of the extracted value of q$, 
defined in eq. ([6| , for the A A system as a function of the number of measurements per configuration. 
The vertical axis is in units of 1/MeV 4 . 



potentials are constructed to reproduce the experimentally-measured NN scattering cross- 
sections. 

The detailed exploration of the behavior of the signal-to-noise ratio in the baryon-baryon 
correlation functions, made possible by the very large number of measurements that have 
been performed, has been exceptionally illuminating. The importance of the "Golden Win- 
dow" of time-slices in which the signal-to-noise ratio is essentially independent of time cannot 
be overstated. This window allows for precise determinations of the energy-splitting between 
interacting baryons and isolated baryons, and in this window, the signal-to-noise ratio does 
not scale with baryon number, making precise measurements in multi-baryon systems feasi- 
ble as discussed in Ref. [23]. 

This calculation is the first part of a thorough analysis of baryon-baryon scattering at this 
pion mass. Calculations in lattice volumes that are both larger and smaller than the present 
lattice volume are underway and will provide measurements of the scattering amplitude 
at two additional momenta. In most channels, this will allow for a determination of the 
scattering parameters (scattering lengths and effective-range parameters) at this pion mass. 
However, it is important to keep in mind that all of these measurements will be at a single 
lattice spacing, b s ~ 0.123 fm. In order to make precise statements, even at this larger pion 
mass, measurements at smaller lattice spacings will be required. 
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FIG. 28: The extracted values of the inverse of the real part of the inverse of the AA scattering 
amplitude versus the extracted value of |k| 2 , defined in eq. (|6j. The measured values correspond 
to different numbers of measurements per configuration. 
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